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Abstract 

Open charm mesons produced in high energy A-A interactions are 
expected to be powerful probes to investigate the medium produced 
in the collision. In this context it is important to measure the pro- 
duction of as many charmed hadrons as possible, such as D , D + , D+ 
and A c , because the measurement of their relative yield can provide 
information on the hadronization mechanism and is necessary to re- 
duce the systematic error on the absolute cross section. The ALICE 
experiment at the LHC is designed to perform such measurements at 
midrapidity down to px below 1 GeV/c, mainly by means of the silicon 
vertex and tracker detector, the time projection chamber and the time 
of flight detector. One of the main channels for the detection of charm 
production in ALICE is the exclusive reconstruction of the D + meson 
through its three charged body decay K~tt + tt + in Pb-Pb (i/s = 5.5 
TeV) and pp (i/i = 14 TeV) collisions. The selection strategies for 
this analysis and the results of a feasibility study on Monte Carlo 
events will be presented together with the perspectives for the study 
of D + quenching and azimuthal anisotropy measurements. 



1 Heavy quarks in nucleus-nucleus collisions 
at high energy 

Heavy-quark production takes place in partonic scatterings during the early 
stages of the nucleus-nucleus collision. The time scale for a cc pair production 
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is ~ h/(2m Q c 2 ) ~ (0.2GeV fm c~ 1 )/(2.4GeV) ~ 0.1 fm/c, which is much 
smaller than the expected lifetime of the Quark Gluon Plasma ~ 10 fm/c. 
Thus, heavy quarks are expected to provide information about the hottest 
initial phase. The measurement of D mesons can be used to extract the charm 
production cross section. The measurement of charm cross section in both pp 
and AA collisions is useful to evaluate the scaling mechanisms which govern 
the charm production from pp to AA collisions. Several nuclear effects can 
break the binary scaling estimated on a geometrical basis with the Glauber 
model p. They are divided into two classes: initial- and final-state effects. 
The former, such as nuclear shadowing, affect heavy-quark production by 
modifying the parton distribution functions in the nucleus. The latter can 
be due to the interactions of the partons in the medium. They affect heavy- 
quark production by modifying the fragmentation function. Energy loss of 
the partons and development of anisotropic flow are among such processes. 

2 Measurements of open charm in the AL- 
ICE experiment 

A benchmark study was performed on the reconstruction of D° mesons in the 
K~n + decay channel [3]. The study of ratios such as D + /D°, D+/D , D+/D+ 
provide information on the hadronization mechanisms. In addition, the re- 
construction of as many charmed states as possible (D°, D + , D+, A C) . . . ) is 
required in order to reduce the systematic errors on the charm cross section. 
The detection of charmed particles is possible in ALICE thanks to the good 
expected performance of the central barrel detectors, namely the Inner Track- 
ing System (ITS), the Time Projection Chamber (TPC) and the Time Of 
Flight detector (TOF). The ALICE barrel is embedded in a magnetic field, 
B=0.5 T, oriented along the beam axis (2-axis). The tracks are curved in the 
transverse plane (rip). The expected Pt resolution is ~ 1%, the primary ver- 
tex resolution is expected to be ~ 10 /im (in case of Pb-Pb) and the expected 
resolution on the minimum rip distance of the prolonged track with respect 
to the primary vertex (impact parameter) is ~ 60 /im for p r ~ 1 GeV/c. 
ALICE provides charged hadron identification from momenta of about 0.3 
GeV/c to about 3 GeV/c. 

2.1 Reconstruction of D + — > K 7r + 7r + decay vertices 

An accurate measurement of secondary vertices is important because the 
weakly decaying charm states have typical decay lengths of the order of hun- 
dreds /xm. The method for the reconstruction of secondary vertices originated 
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by 3-prong decays is based on the straight line approximation of the tracks, 
i.e. each track which is curved by the magnetic field is approximated by a 
straight line in the vicinity of the primary vertejJl]. The algorithm finds the 
point of minimum distance between three straight lines by minimizing the 
quantity D 2 = S j=li2i3 (x ^ o)2 + i«£ + (*-*>)" ; where (x ,y ,z ) are the 

xi yi zi 

coordinates of the secondary vertex and a x , a y , a z are the errors on the track 
parameters. At high pr of the D + , its decay products get more and more 
collinear with the D meson p?. Hence, along this direction a worsening of the 
resolution of the secondary vertex finder should be observed. This behaviour 
is illustrated in the left panel of Fig. HJ where the RMS is calculated in a 
reference system in the bending plane in which x is the direction of the pt of 
the D + and y is the axis orthogonal to x. Along the x coordinate the RMS 
increases at high px, while in the y coordinate the RMS goes down to ~ 40 
/iin. 

3 Selection strategy for the D + — > K~tt + tt + chan- 
nel in Pb-Pb collisions. 

The expected number of D + in a central Pb-Pb event is ~ 4 in \y\ < 1 [3]- 
We generated separately signal and background events. The results pre- 
sented in this paper were obtained with a sample of 20 000 central HI- 
JING (dA^ c /j/d|/=6000) events and ~ 2 x 10 6 signal triplets generated with 
PYTHIA. The signal events have approximately the same particle multiplic- 
ity as the central HIJING events. A sample of ~ 5.4 x 10 6 minimum bias pp 
events was generated with PYTHIA; due to the smaller CPU time needed to 
simulate and reconstruct pp events compared to Pb-Pb, a separate generation 
of signal and background events was not necessary. The cuts presented in 
the following are optimized in order to maximize the significance S/ V 'S + B 
in the invariant mass range \M inv — Mjj±\ < la = 12 MeV/c 2 . The strategy 
is based on a series of selection steps. The first selection is done on the sin- 
gle tracks before they are combined into triplets and passed to the secondary 
vertex finding algorithm. The goal is to reduce the number af all the possible 
combinations of tracks into triplets in a central Pb-Pb collision. This number 
is ~ 10 10 without any initial cut and it must be cut before using the vertexer. 
The best combination of cuts which preserves the px spectrum of D ± mesons 
down to ~ 1 GeV/c is defined by prcut,K = 0.7 GeV/c, pr C ut,-K = 0.5 GeV/c 
and d 0tCUt = 95 jttm (do is the distance of closest approach of the track to 
the primary vertex in the bending plane); with such a choice, the number 

1 We verified that the error coming from the straight line approximation is negligible. 
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Figure 1: Left panel: RMS of the residual distribution for the x,y and z 
coordinates of the secondary vertex as a function of the pt of the D + . Right 
panel: ratio of triplets (both signal and background) selected by the cut 
a < cr cu t, for D ± transverse momentum 2 < < 3 GeV/c. 



of combinatorial background triplets per event is reduced to ~ 10 . The 
second selection takes place before the triplets are passed to the vertexer 
and it is based on a cut on the distance 5 between the primary vertex and 
the secondary vertex of two opposite charged tracks. The chosen value of 
5 cut is ~ 700 /im. The triplets satisfying specific requirements on combina- 
tions of track impact parameters are then selected and the secondary vertex 
finding is applied. A selection based on the track dispersion is performed 
(a 2 = £1=1,2,3(2^ - xo) 2 + (Vi - Vo) 2 + {zi - ^o)€l)- The right-hand plot in 
Fig. [I] reports the fraction of accepted triplets (both signal and background) 
with a < a cut , for D 1 * 1 transverse momentum 2 < p T < 3 GeV/c. The value 
a cut has been tuned in various intervals of D ± transverse momentum; it is of 
the order of 200 fim. The candidate triplets undergo additional selections. 
The cut variables which have been considered at this level of the selection 
strategy are the following: 

1. Distance between the primary and the secondary vertices, d. The signal 
triplets are characterized by larger values of d, according to the fact 
that they originate from displaced decay vertices 

2. Maximum transverse momentum among the three tracks belonging to a 
triplet, indicated as Max{p TK ,p Tn i,p T7T2 }- This cut accounts for the 
slightly harder p? of the decay products (pions and kaons) of D mesons 
compared to all the pions and kaons produced in a Pb-Pb event. 

2 No weight with the errors on the track parameters was used; in this way a has the 
meaning of a standard deviation. 
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3. cos6 point , where 9 point is the pointing angle, defined by the direction of 
the reconstructed D 1 * 1 momentum on the bending plane and the line 
connecting the primary and the secondary vertices. If the found vertex 
really corresponds to a D 1 * 1 decay vertex, then O po i n t ~ and cos 9 po i nt ~ 
1. 

4. Sum of the squares of the impact parameters of the three tracks with re- 
spect to the primary vertex, s 2 = 1^=1,2,3 <^o % = <^o k + <^o ^ + <^o n 2 ■ This 
variable is expected to be larger for signal triplets which are displaced 
from the primary vertex. 

A method which simultaneously tunes all the four cut variables described 
above was developed. This method is based on the creation of multi-dimensional 
matrices for both signal and background events. Each cell ijkl of the multi- 
dimensional matrices of signal and background contains the number of triplets 
which pass the selection criteria defined by the set of cuts: Cijki — {[d > 
d cu u] AND [p M > PMcutj] AND [cos 9 point > cos 9 pointiCUt k ] AND [Ed^ > 
(Sc?q ^cutj}}- This procedure is applied in several intervals of D 1 * 1 transverse 
momentum p?. Three analyses were performed in Pb-Pb, namely with ideal 
PID, with real PID and without any PID. The D 1 * 1 reconstruction was also 
studied in pp collisions assuming no particle identification. 



4 Results 

An example of a result of the optimization of the cuts with the four-dimensional 
matrix procedure for 2 < < 3 GeV/c is reported in Fig. [2, where the 
matrix of the significance is projected on the plane (cos^o^^, d cut ), for 
and PM,cut = 0.5 GeV/c. The significance S/y/S + B 
maximized in the D transverse momentum bins < pt < 2 GeV/c, 2 < 
Pt < 3 GeV/c, 3 < pt < 5 GeV/c, pr > 5 GeV/c is reported in Fig. [H] for 
both Pb-Pb collisions (for the three cases of ideal PID, real PID and without 
PID) and pp collisions (without PID). The obtained results were used to eval- 
uate the statistical uncertainities related to the measurement of the nuclear 
modification factoiQ and of the v 2 coefficient [3] , which quantifies the amount 
of azimuthal anisotropy in the pt distribution. The errors are evaluated on 
a statistics of D 1 * 1 mesons expected for a one year run. Raa is calculated [1] 
as a function of pt in the case of no particle identification for both Pb-Pb 
and pp collisions and is shown in the left-hand panel of Fig. H] for different 

3 This cut is not useful in this particular pt bin. 

4 The nuclear modification factor is defined as: BR a = . V?d Aa £i — 

AA yicld-y ' xN con 
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Figure 2: Significance (for 10 7 Pb-Pb central events) as a function of the cut 
variables d cut and cos8 point ,cut, for p M ,cut = 0.5 GeV/c and (Edj^)^ = /zm 2 . 
The transverse momentum of the triplets is 2 < px < 3 GeV/c. The real 
PID case is considered. 
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Figure 3: Left panel: significance S/y/S + B (normalized to 10 7 central Pb- 
Pb events) as a function of the D 1 * 1 transverse momentum for ideal PID 
(black circles), real PID (red triangles up) and no PID (blue triangles down). 
Right panel: significance (normalized to 10 9 minimum bias pp collisions) as 
a function of the D 1 * 1 transverse momentum for pp collisions without PID. 
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values of the transport coefficient. The performance for the measurement of 
the i>2 coefficient was estimated with a fast simulation consisting in the ran- 
dom generation of D mesons distributed with an azimuthal anisotropy. The 
number of D ± depends on pt and on centrality. The generated D azimuthal 
angle is smeared with its experimental resolution. Each of the generated D 1 * 1 
is superimposed to an event made of tracks generated with a vf coefficient 
(different from the D 1 * 1 meson v®)- These particles are used to reconstruct 
the event plane. The results of v^ipr) for D mesons are reported in Fig. H] 
(right-hand panel); the values of of D 1 * 1 mesons used as input of the sim- 
ulation are taken from [6J. The predictions from [6] are calculated in the 
framework of a coalescence model in case of charm flow (solid line) and in 
case of flow of only light quarks (dashed dotted line). 



Figure 4: Left panel: The Raa of D mesons as a function of px, calculated 
with the parameter q = 25 GeV 2 /fm (dark grey curve). The statistical and 
systematic errors are described in the legend. The results are compared to 
the scenario without quenching of D 1 * 1 (q = GeV 2 /fm, light grey curve) 
and with q = 100 GeV 2 /fm including the dead cone effect (dashed-dotted 
light grey line). Right panel: D meson v^ipr) as a function of pt in the 
centrality class 6 < b < 9 fm for a choice of the transport model parameter 
Up = 10 mb. 

In conclusion, the reconstruction of the channel D + — ► K~7r + 7r + seems 
to be feasible with a good significance both in Pb-Pb and in pp collisions. 
In the likely case in which the multiplicity in central Pb-Pb events turned 
out to be lower, the significance would of course improve. The analysis on 
Raa is feasible in one year of data taking. The statistical error bars for the 
measurement of function of the transverse momentum in one year 

of data taking are a bit large. The use of a semi-peripheral trigger and the 
use of the D meson statistics collected in several years would improve the 
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situation. 
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